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We present the final combination of CDF and DO measurements of cross sections for single-top- 
quark production in proton-antiproton collisions at a center-of-mass energy of 1.96 TeV. The data 
correspond to total integrated luminosities of up to 9.7 fb~^ per experiment. The t-channel cross 
section is measured to be at = pb. We also present the combinations of the two-dimensional 

measurements of the s- vs. t-channel cross section. In addition, we give the combination of the s + t 
channel cross section measurement resulting in as+t = S.SOjlQ ^g pb, without assuming the standard- 
model value for the ratio as fat- The resulting value of the magnitude of the top-to-bottom quark 
coupling is \Vtb\ = 1.0210;°®, corresponding to \Vtb\ > 0.92 at the 95% C.L. 

PACS numbers: 14.65.Ha; 12.15.Ji; 13.85.Qk; 12.15.Hh 


The top quark is the heaviest elementary particle of 
the standard model (SM). Detailed studies of top-quark 
production and decay provide stringent tests of strong 
and electroweak interactions, as well as sensitivity to ex¬ 
tensions of the SM [l[. At the Fermilab Tevatron collider, 
protons (p) and antiprotons (p) collided at a center-of- 
mass energy of y/s = 1.96 TeV. Top quarks were pro- 
ducedpredominantly in pairs {tt) via the strong interac¬ 
tion [2|. They were also produced singly via the elec¬ 
troweak interaction. The cross section for single-top- 
quark production depends on the square of the magni¬ 
tude of the quark-mixing Cabibbo-Kobayashi-Maskawa 
(CKM) matrix Q element Vtb, and consequently is sensi¬ 
tive to contributions from a fourth family of quarks in, 
as well as other new phenomena @ , which would lead to 
a measured strength of the Wtb coupling \Vth\ different 
from the SM prediction. Non-SM phenomena could also 
change the relative fraction of events produced in the 


various channels that contribute to the total single-top- 
quark production cross section. 

In pp scattering, single-top-quark production proceeds 
in the t-channel via the exchange of a spacelike virtual 
W boson between a light quark and a bottom quark [3- 
1^. Single top quarks are also produced in the s chan¬ 
nel via the decay of a timelike virtual W boson pro¬ 
duced by quark-antiquark annihilation, which produces 
a top quark and a bottom quark 0 , or in association 
with a W boson {Wt) 0- The predicted SM cross sec¬ 
tion for the t-channel process at is 2.10 ± 0.13 pb |9|, 
while the s-channel cross section Ug is 1.05±0.06 pb [12| . 
both calculated at next-to-leading-order (NLO) in quan¬ 
tum chromodynamics (QCD), including next-to-next-to- 
leading logarithmic (NNLL) corrections. A top-quark 
mass of 172.5 GeV was chosen, which is consistent with 
the current world average value 0. The cross section 
for Wt production awt is negligibly small at the Teva- 
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tron and therefore is not considered in the combination 
described in this Letter. Since the magnitude of the Wtb 
coupling is much larger than that of Wtd or of Wts 0, 
each top quark decays almost exclusively to a W boson 
and a b quark. 

Observation of si ngle-top -quark production was re¬ 
ported by the CDF [iMl and DO [H 0 Collabora¬ 
tions not differentiating between the s and the t chan¬ 
nels (hereinafter s + t channel). The CDF Collabora¬ 
tion subsequently measured a single-top-quark produc¬ 
tion cross section for the sum of the s, t, and Wt chan¬ 
nels of (Js+t+wt = 3.04+q;| 3 pb using data correspond¬ 
ing to 7.5 fb“^ of integrated luminosity [0 and for the 
sum of the s and t channels of ag+t = 3.02^;4g pb us¬ 
ing up to 9.5 fb“^ of integrated luminosity [2l|. The DO 
Collaboration obtained cTs+t = pb using data 

corresponding to 9.7 fb“^ of integrated luminosity |22|. 

The cross sections for individual production modes 
were also measured separately. The DO Collaboration ob¬ 
served the t-channel process 0 and measured its cross 
section to be at = 3.07to;4g pb using data correspond¬ 
ing to 9.7 fb“^ of integrated luminosity [0. The CDF 
Collaboration measured at+wt = 1-661q 47 using data 
corresponding to 7.5 fb“^ of integrated luminosity (0 
and at = l-65lQ;3g P^ using up to 9.5 fb“^ [ilj of in¬ 
tegrated luminosity. The difference between the results 
for at is about 2 standard deviations (s. d.). Further¬ 
more, both the CDF and DO Collaborations reported 
evidence for s-channel production 22, 24, 0 and com¬ 
bined their results to observe the s-channel process with 


= 1.29 


-1-0.26 

-0.24 


pb 


At the CERN LHC proton-proton {pp) collider, t- 
channel production was observed by the ATLAS and 


CMS Collaborations [27H30j. Furthermore, ATLAS has 
found evidence for Wt associated production 3l|, fol¬ 
lowed recently by an observation at the CMS experi¬ 
ment [s^. All measurements are in agreement with SM 
predictions i0. 

In this Letter, we report final combinations of single- 
top-quark cross section measurements from analyses per¬ 
formed by the CDF 211 and DO [0 Collaborations using 
up to 9.7 fb“^ of integrated luminosity per experiment. 
In particular, we present a combined t-channel cross sec¬ 
tion, a combined two-dimensional measurement of the 
s- vs t-channel cross sections, and a combination of the 
(s -I- t)-channel cross sections. The combination is ob¬ 
tained by collecting the inputs from both experiments 
and reperforming the statistical analysis. This approach 
allows for a tighter constraint on the systematic uncer¬ 
tainties that are common to both experiments, leading 
to a higher precision than that achievable from aver¬ 
aging the individual results. Here, we do not include 
the combination of the s-channel cross-section measure¬ 
ments, which was reported in Ref. [0. We also measure 
the magnitude of the CKM matrix element Vtb with no 
assumptions on the number of quark flavors. 


The CDF and DO detectors are large solenoidal 
magnetic spectrometers surrounded by projective-tower- 
geometry calorimeters and muon detectors [3i,|3i|. The 
data were selected using a logical OR of many online se¬ 
lection requirements that preserve high signal efficiency 
for offline analysis. Both collaborations analyze events 
with a lepton (£ = e or p) plus jets and an imbalance 
in the total event transverse energy reconstructed as 
the negative vector sum of all significant transverse en¬ 
ergies in the calorimeter cells and the muon transverse 
momenta subtracting the calorimeter energy deposition 
due to muons (£-|-jets). This topology is consistent with 
single-top-quark decays in which the decay W boson sub¬ 
sequently decays to ^v 0 [ 23 . Events were selected 
that contain only one isolated lepton i with large trans¬ 
verse momentum large I^lT^ and two or three clus¬ 
ters of energy in the calorimeters (jets) with large pt- 
One or two of these jets were required to be identified as 
emerging from the hadronization of a 6 quark (5-tagged 
jets). Multivariate techniques were used to discriminate 
&-quark jets from light-quark and gluon jets [3l,|33. Ad¬ 
ditional selection criteria were applied to exclude kine¬ 
matic regions that were difficult to model and to min¬ 
imize the background of multiple jets from QCD pro¬ 
duction (QCD multijet) in which one jet was misrecon- 
structed as a lepton and spurious arose from mismea- 
surements 0.[^. 

The other final-state topology, analyzed by the CDF 
Collaboration, involves jets, and no reconstructed 
isolated charged leptons (fi^T+jets) [ilj. In the CDF 
.^'r+jets analysis, overlap with the .^-l-jets sample was 
avoided by vetoing events with identified leptons. Large 
was required, and events with either two or three re¬ 
constructed jets were accepted. This additional sample 
increased the acceptance for signal events by including 
those in which the W-boson decay produced a lepton that 
is either not reconstructed or not isolated, or a r lepton 
that decayed into hadrons and a neutrino, which were re¬ 
constructed as a third jet. After the basic event selection, 
QCD multijet events dominate the £?T+jets event sam¬ 
ple. To reduce this background, a selection based on an 
artificial neural network was optimized to preferentially 
select signal-like events 211. 

Events passing the £-|-jets and£^T+jets selections were 
separated into independent channels based on the num¬ 
ber of reconstructed jets as well as on the number and 
quality of 5-tagged jets. Each of the channels has a dif¬ 
ferent background composition and signal-to-background 
ratio, and analyzing them separately enhances the sen- 
sitivity to single-top-quark production by approximately 

10% [illiil. 

Several differences in the properties of s- and t-channel 
events were used to distinguish them from one another. 
Events originating from f-channel production typically 
contain one light-flavor jet at large pseudorapidity mag¬ 
nitude \r]\, which is useful for separating them from events 
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associated with s-channel production and other SM back¬ 
ground processes. Events from the s-channel process are 
more likely to yield two b jets within the central region 
of the detector. 

Both collaborations used Monte Carlo (MC) event gen¬ 
erators to simulate kinematic properties of signal and 
background events, except for multijet production, which 
was modeled with data using matrix methods [371, [SSj. 
Using the powheg generator, CDF modeled single- 
top-quark signal events at NLO accuracy in the strong 
coupling strength as- This is different from DO where the 
SINGLETOP event generator was used, based on NLO 
QCD COMPHEP calculations that match the kinematic 
features predicted by other NLO calculations ^ l42 |. 
Spin information in the decays of the top quark and the 
W boson is preserved in both powheg and SINGLETOP. 

Kinematic properties of background events from pro¬ 
cesses in which a W or Z boson is produced in associ¬ 
ation with jets (fU-l-jets or Z-l-jets) were simulated us¬ 
ing the ALPGEN MC generator for the calculation of 
tree-level matrix elements interfaced to PYTHIA ^ for 
parton showering and hadronization and using the MLM 
matrix-element parton-shower matching scheme ^ . Di¬ 
boson contributions (WW. WZ, and ZZ) were mod¬ 
eled using PYTHIA The ti process was modeled 

using PYTHIA at CDF and ALPGEN at DO. The mass 
of the top quark in simulated events was set to rrit = 
172.5 GeV. Higgs-boson processes were modeled using 
simulated events generated with pythia for a Higgs bo¬ 
son mass of rriH = 125 GeV 3“ 4^. In all of the above 
cases, PYTHIA was used to model proton remnants and 
to simulate the hadronization of all generated partons. 
The presence of additional pp interactions was modeled 
by overlaying events selected from random beam cross¬ 
ings matching the instantaneous luminosity profile in the 
data. All MC events were processed through GEANT- 
based detector simulations |49l |. and were reconstructed 
using the same computer programs as used for data. 

Data were used to normalize IV-boson production asso¬ 
ciated with both light- and heavy-flavor jet contributions 
in samples enriched in VF-l- iets processes, which have neg¬ 
ligible signal content [iTI, l22|, [23- All other simulated 
background samples were normalized to their theoreti¬ 
cal cross sections, i.e., ti at next-to-next-to-leading or¬ 
der QCD [H^, Z+iets and diboson production at NLO 
QCD [^, and Higgs-boson production including all rele¬ 
vant higher-order QCD and electroweak corrections . 
For the measurement of at, the s-channel single-top- 
quark production sample was considered as background 
and normalized to the NLO QCD cross section combined 
with NNLL resummations 121. 

Multivariate discriminants were optimized to separate 
signal events from large background contributions. To 
combine the results from the two experiments, we use 
the s- and t-channel discriminants from the CDF 
and DO [1^ single-top-quark measurements. We per¬ 


form a likelihood fit to the binned distribution of the 
final discriminants. We combine the varions channels of 
the different analyses from each experiment by taking 
the product of their likelihoods and simultaneously vary¬ 
ing the correlated uncertainties and by comparing data 
to the predictions for each contributing signal and back¬ 
ground process. Using a Bayesian statistical analysis (53| . 
we then derive combined Tevatron cross section measure¬ 
ments, taking the prior density for the signal cross sec¬ 
tions to be uniform for non-negative cross sections. 

For the sources of uncertainties we follow Ref. [26| . 
We consider the following systematic uncertainties: the 
integrated luminosity from detector-specific sources and 
from the inelastic and diffractive cross sections. We also 
consider systematic uncertainties on the signal model¬ 
ing, the simulation of background, data-based methods 
to estimate background, detector modeling, 5-jet tagging, 
and the measurement of the jet-energy scale. Table I of 
Ref. [ 2 ^ summarizes the categories that contribute to the 
uncertainties on the shape of the output of the multivari¬ 
ate discriminants distributions and the range of uncer¬ 
tainties applied to the predicted normalizations for sig¬ 
nal and background contributions. Reference gives 
the sources of systematic uncertainty common to mea¬ 
surements of both collaborations that are assumed to be 
fully correlated, and lists uncertainties that are assumed 
to be uncorrelated. The dependence of the results on 
these correlation assumptions is negligible. 

A two-dimensional (2D) posterior-probability density 
is constructed as a function of as and at in analogy to 
the one-dimensional (ID) posterior probability described 
in Ref. [i^. The measured cross section is quoted as 
the value at the position of the maximum, and the 68% 
probability contour defines the measurement uncertainty. 

Figure[T]shows the distribution of the mean values from 
the discriminants sorted by the s-channel minus t-channel 
expected signal contributions divided by the background 
expectation, (s — t)/b. An entry in the histogram corre¬ 
sponds to a collection of bins with similar ratio (s — i)/b. 
The value on the horizontal axis is given by the mean 
discriminant for those bins. The vertical axis gives the 
number of events in those bins. We show the data, the 
SM predictions for the s- and t-channel processes, and 
the predicted backgrounds separated by source. The dis¬ 
tribution for large negative values is dominated by the 
content of the bins that show a higher t-channel contri¬ 
bution, while large positive values are dominated by the 
content of the bins with a higher s-channel contribution. 
The abscissa extends to larger negative values since we 
expect more t-channel events than s-channel events and 
the separation from background is better for t-channel 
events than for s-channel events. The region correspond¬ 
ing to discriminant values near zero is dominated by the 
background. 

FigurelUpresents the resulting 2D posterior probability 
distribution as a function of at and as- The value and 










Discriminant output [(s-channel - t-channel)/background] 


FIG. 1: Distribution of the mean discriminants for bins 
with similar ratios of (s-channel—f-channel) signals divided by 
background yields. The data, predicted SM s- and t-channel 
yields, and expected background are displayed. The total 
expected background (black solid line) is shown with its un¬ 
certainty (gray shaded band). A nonlinear scale is used on 
the abscissa to better display the range of the discriminant 
output values. 


uncertainty in the individual cross sections are derived 
through the ID posterior probability functions obtained 
by integrating the 2D posterior probability over the other 
variable. The most probable value of at is 2.251 q3 ® pb. 
The measurement of is performed without making 
assumptions on the ratio of as I at by forming a 2D poste¬ 
rior probability density distribution of ag+t versus at and 
then integrating over all possible values of at to extract 
the ID estimate of ag+t- The combined cross section is 
CTs+t = 3.30 ('(q |q pb. The total expected uncertainty on 
CTs+t is 13%, the expected uncertainty without consid¬ 
ering systematic uncertainties is 8%, and the expected 
systematic uncertainty is 10%. The systematic uncer¬ 
tainty from the limited precision of top-quark mass mea¬ 
surements is negligible [l7|,l23. Figure [2] also shows the 
expectation from several beyond the SM (BSM) models. 
Figure [3] shows the individual [2l|,[22j and combined (this 
Letter) measurements of the t- and (s -I- t)-channel cross 
sections including previous measurements of the individ¬ 
ual [i^ [ 2 ^ and combined [ 2 ^ s-chaimel cross sections. 
All measurements are consistent with SM predictions. 

The SM single-top-quark production cross section is 
directly sensitive to the square of the CKM matrix ele¬ 
ment Vtb (ini, thus providing a measurement of \Vtb\ 
without any assumption on the number of quark families 
or the unitarity of the CKM matrix [s^. We extract \ Vtb\ 
assuming that top quarks decay exclusively to Wb hnal 
states. 

We start with the multivariate discriminants for the 
s and t channels for each experiment and form a 
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FIG. 2: Two-dimensional posterior probability as a function 
of at and as with one s. d. (68% G.L.), two s.d. (95% 
C.L.), and three s. d. (99.7% C.L.) probability contours 
for the combination of the CDF and DO analysis channels 
compared with the NLO+NNLL theoretical prediction of the 
SM d, [T 3 . Several BSM predictions are shown, a model 
with four quark families with top-to-strange quark coupling 
|Fts| = 0.2 d], a top-flavor model with new heavy bosons with 
mass rrix = 1 TeV [3, a model of charged top pions with mass 
m^± — 250 GeV [y|, and a model with flavor-changing neu¬ 
tral currents with a 0.036 coupling Ku/A between up quark, 
top quark, and gluon Bin. 


Bayesian posterior probability density for \Vtb\'^ assum¬ 
ing a uniform-prior probability distribution in the region 
[0, 00 ] corresponding to a uniform prior density of the 
signal cross section. Additionally, the uncertainties on 
the SM predictions for the s- and t-channel cross sec¬ 
tions d,[l3 are considered. The resulting posterior prob¬ 
ability distribution for \Vtb\^ is presented in Fig. |T| We 
obtain \Vtb\ = 1.02t°;"t If we restrict the prior to the 
SM region [0,1], we extract a limit of \Vtb\ > 0.92 at the 
95% C.L. 

In summary, using pp collision samples corresponding 
to an integrated luminosity of up to 9.7 fb“^ per experi¬ 
ment, we report the final combination of single-top-quark 
production cross sections from CDF and DO measure¬ 
ments assuming mt = 172.5 GeV. The cross section for 
t-channel production is found to be 

at = 2.25t°o fi pb. 

Without assuming the SM value for the relative s- and 
t-channel contributions, the total single-top-quark pro¬ 
duction cross section is 

ag+t = 3.30l°ig pb. 
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Tevatron Run II single top quark summary 
Measurement Cross section [pb] 
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FIG. 3: Measured single-top-quark production cross sections 
from the CDF and DO Collaborations in different produc¬ 
tion channels and the Tevatron combinations of these analyses 
compared with the NLO-I-NNLL theoretical prediction [9l.[l2|. 
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FIG. 4: Posterior probability distribution as a function of 
\Vtb\^ for the combination of GDF and DO analysis channels. 
The arrow indicates the allowed values of |Vi6|^ corresponding 
to the limit of \Vtb\ > 0.92 at the 95% G.L. 


Together with the combined s-channel cross section , 
this completes single-top-quark cross-section measure¬ 
ments accessible at the Tevatron. All measurements 
are consistent with SM predictions UM- Finally, we 
extract a direct limit on the CKM matrix element of 
\Vtb\ > 0.92 at the 95% C.L. As a result, there is no 
indication of sources of new physics beyond the SM in 
the measured strength of the Wtb coupling. 
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